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Abstract 

r*j ' We consider the asymptotics of the correlation functions of the characteristic 

Q-e polynomials of the hermitian Wigner matrices H n = n~ x l 2 W n . We show that for 

the correlation function of any even order the asymptotic coincides with this for the 
GUE up to a factor, depending only on the forth moment of the common probability 
law Q of entries QWj k , ^RWj k , i.e. that the higher moments of Q do not contribute 
to the above limit. 

7— I ■ 

vo : 1 Introduction 

CO 

Characteristic polynomials of random matrices have been actively studied in the last years. 
. The interest was initially stimulated by the similarity between the asymptotic behavior of 
the moments of characteristic polynomials of a random matrix from the Circular Unitary 
Ensemble and the moments of the Riemann (^-function along its critical line (see jH]). But 
with the emerging connections to the quantum chaos, integrable systems, combinatorics, 
representation theory and others, it has become apparent that the characteristic polyno- 
mials of random matrices are also of independent interest. This motivate the asymptotic 
study of the moments of characteristic polynomials for other random matrix ensembles 
(see e.g. [10], [3]). 

In this paper we consider the hermitian Wigner Ensembles with symmetric entries 
distribution, i.e. hermitian n x n random matrices 

H n = n- l ' 2 W n (1.1) 

with independent (modulo symmetry) and identically distributed entries 3lWj tk and ^sWj )k 
such that 

E{W jk } = E{(W jk ) 2 } = 0, E{|W>| 2 } = 1, n v 

E{^ 2l+1 W jk } = E{S 2l+1 W jk } = 0, j, k = 1, .., n, I EN. { ' 

Denote by , . . . , An™"* the eigenvalues of random matrix and define their Normalized 
Counting Measure (NCM) as 

» 



O 



N n (A) = (l{Af G A, j = l,..,n}/n, N n (R) = 1, (1.3) 



1 



where A is an arbitrary interval of the real axis. The global regime of the random matrix 
theory, centered around the weak convergence of the Normalized Counting Measure of 
eigenvalues, is well-studied for many ensembles. It is shown that N n converges weakly to 
a non-random limiting measure N known as the Integrated Density of States (IDS). The 
IDS is normalized to unity and is absolutely continuous in many cases 

N(R) = 1, N(A) = [ p(X)dX. (1.4) 



The non-negative function p in (11.41) is called the limiting density of states of the ensemble. 
In the case of Wigner hermitian ensemble it is well-known (see, e.g.,[H]) that 



p(A) = Psc (X) = —VA^V. (1.5) 
The mixed moments or the correlation functions of characteristic polynomials are 

„ 2m 

^ 2m (A) = J UdetiX, - H)P n {dH n ), (1.6) 



tin 



where H n is the space of hermitian n x n matrices, 

n 

dH n = ]JdH JJ J] 1iV/,,a// ; , (1.7) 

j=l ^<j<k<n 

is the standard Lebesgues measure on "H n , P n (d H n ) is probability law of the nxn random 
matrix H n , and A = {Ay}^ are real or complex parameters that may depend on n. 
We are interested in the asymptotic behavior of (11.61) for matrices (II. ip as n — > oo for 

Aj = A H 7T~T' J = 1 >--> 2m > 

^P SC (A ) 

where A G (—2,2), p sc is defined in (I1.5P and £ = are real number varying in a 

compact set K C R. 

In the case of hermitian matrix model, i.e. the matrices with 

P n (dH n ) = Z- 1 e- ntTV ^dH n , 

where V is a potential function, the asymptotic behavior of (II. 6p is known. Using the 
method of orthogonal polynomials, it was shown (see [13], [2]) that 

' F 2m f A Q + e/(np(A )) 



2m 



mnV(\ )+a v (\ ) E -„ . ,„ 

where A = (A , . . . , A ) G M 2m , 

av{X) = 2p\X) ' 
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p is a density of (11.41) . A is such that p(A ) > and A(xi, . . . ,x m ) is the Vandermonde 
determinants of X\, . . . , x m . 

Unfortunately, the method of orthogonal polynomials can not be applied to the general 
case of hermitian Wigner Ensembles. Thus, to find the asymptotic behavior of (II. 6ft other 
methods should be used. In f7j Gotze and Kosters use the exponential generating function 
to study this behavior for the second moment, i.e. for the case m = 1 in (II .6p . In this 
case it was shown for matrices (II. ip that 

-4rr^2 (Ac + £i/(np sc (\ )), A + 6/MAo))) 
np{A ) 

= 2neMn(Xl ~ 2)/2 + a(A )(6 + 6) + I «* 1 } ^^St^ ( 1 + 

where 

a ( A ) = o — mv k 4 = /^4-3/4, (1.8) 
2p sc (A) 

and /i4 is the forth moment of the common probability law Q of entries QWjk, 9ftW}&. 

In this paper we consider the general case m > 1 of ( 11. 6ft for the random matrices 
(HI}. Set 

£ (n) (0 = — TTT-^ 2 f A ° 7TT' A ° TTtI 

np(A ) V ^c(Ao) ^p sc (A )y 

= 27rexp {^(A 2 - 2) + 2a(A )£ + ^} (1 + o(l)). 

The main result of the paper is 

Theorem 1. Lei £/ie entries ^sWjk, 3tWjk of matrices U.l\) has the symmetric probability 
distribution with 4m finite moments. Then we have for m > 1 

hm — F 2m C A + £/(np sc (\ )) 

(^ sc (A )) m2 n v^te) 
i=i 

= exp{m(m - l)/t 4 (A§ - 2) 2 /2} f sin(7rfe - £ m+J -)) ) ? 

^M)A(a,..,gA(u,.,6 w ) e I vr^-w J, 

where F 2m and p sc (A) are defined in U.6\) and U.5\) . A = (A , . . . , A ) G M 2m ; A G 
(—2,2), £ = {£j} 2 ™ 1; and a(A) and k 4 are defined in U.8\) . 

The theorem shows that the above limit for the mixed moments of characteristic 
polynomials for random matrices (11.11) coincide with those for the GUE up to a factor, 
depending only on the forth moment of the common probability law Q of entries QWjk, 
dtWjk, i.e. that the higher moments of Q do not contribute to the above limit. This is a 
manifestation of universality of the limit, that can be composed with universality of the 
local bulk regime for Wigner matrices (see [5]). 

The paper is organized as follows. In Section 2 we obtain a convenient integral rep- 
resentation for F 2m in the case of symmetric probability distribution of entries with 
Am finite moments by using the integration over the Grassmann variables and Harish 
Chandra/Itzykson-Zuber formula for integrals over the unitary group. In Section 3 we 
prove Theorem [1] by applying the steepest descent method to the integral representation. 

We denote by C, Ci, etc. and c, ci, etc. various n- independent constants below, which 
can be different in different formulas. Integrals without limits denote the integrals over 
whole real axis. 
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2 The integral representation. 



In this section we obtain the integral representation for the correlation functions F 2m of 
( II. 6p of characteristic polynomials. To this end we use the integration over the Grassmann 
variables. The integration was introduced by Berezin and widely used in the physics 
literature (see [1] and [1]). For the reader convenience we give an outline of this technique 
here. 

Let us consider the two sets of formal variables {ipj}™ =1 , {i^j}]=i, which satisfy the 
following anticommutation conditions 

ipjipk + ^k^j = ^j?Pk + i/Jkipj = ipjipk + ^k^j = 0) Ji k = l,..,n. 
In particular, for k = j we obtain 

These two sets of variables {■?/>,,• }™ =1 and { , J }" =1 generate the Grassmann algebra A. 
Taking into account that ip? = 0, we have that all elements of A are polynomials of 
{■?/>,• }™ =1 and {'0 J }" =1 . We can also define functions of Grassmann variables. Let x be an 
element of A. For any analytical function / by f(x) we mean the element of A obtained 
by substituting x m the Taylor series of / near zero. Since x is a polynomial of {ipj}^ = i, 
{"0j}^=i? there exists such I that x l = 0) an d hence the series terminates after a finite 
number of terms and so f(x) £ A. 

Following Berezin [TJ, we define the operation of integration with respect to the anti- 
commuting variables in a formally way: 



J dtpj = J dipj — 0, J tpjd ipj = J ipjd ipj 



1. 



This definition can be extend on the general element of A by the linearity. A multiple 
integral is defined to be repeated integral. The "differentials" dipj and dijj k anticommute 
with each other and with the variables ipj and ip k . 
Therefore, if 

m 

f{Xi, ■ ■ ■ , Xm) = «o + ^2 a hXh + a hhXhXh + ■■■ + a 1:2 ,...,mXi ■■■Xm, 

ii=i ji<h 

then 

f(Xi,---iXm)dXm---dxi = ai,2,...,v< 



Let A be an ordinary hermitian matrix. The following Gaussian integral is well-known 



/ 



i \ - , i -A- d ffizjd ^sZj 1 
exp - J2 A hkZ fZk II ^ = T^Ta- ^ 



it detA 

j,k=l ) 3=1 



One of the important formulas of the Grassmann variables theory is an analog of formula 
(12. ip for Grassmann algebra (see pQ): 



/( n ~\ n 

exp <j ^ Aj^ipji/Jk > Y^dtpjdipj = detA, 

lj,k=l J j=l 



(2.2) 
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where {^j}™ =1 and {'0,}?=i are the Grassmann variables. Besides, we have 

/q ( n ~\ n 

IJ^p^p ex P \ A i^j^k > Yl dipjdtpj = detA ht .. Usit .. tSq , (2.3) 

p=l \j,k=l ) j=l 

where Ai lt __ t i q;sij __ jSq is a (n — q) x (n — g) minor of the matrix A without rows l\, .., l q and 
columns si, .., s q . 

2.1 Asymptotic integral representation for F2. 

In this subsection we obtain the asymptotic integral representation of (I1.6P for m = 1. The 
corresponding asymptotic formula was obtained in [7] by using the exponential generating 
function. We give here a detailed proof based on the Grassmann integration to show the 
basic ingredients of our technique. The technique will be elaborated in the next subsection 
to obtain the asymptotic integral representation of (11. 6p for m > 1. 
Set 



D2 = HJD(n)&), (2.4) 



1=1 



where D^^) is defined in (11.91) . Note also that 



^W(() = e »(AoKW2 >( 4 A o) / ( t -a /2)"e-t(wV2) 2 dt(l + o(l)) (2.5) 



as n 00. 

Using (12. 2p . we obtain from ( 11. 6ft 



2 n 2 n 



L> 2 - X F 2 (A) = D^E ^ / e<— 1 JJJJdV^d^ 

I"' r=l 9=1 

f2 n _ /- 2 

/" E E V'psV'i.s 1; /— — /n fi^ 

ye- 1 «pj-x;-^x)(^«+^ii) ( 2 - 6 ) 



j<k v !=1 

zSrau 2 n w 2 I 2 n 

- e — 7^ yicv'ii^w - V'wV'ij) - 5^ -r^yz^ji^i \ n n d^ qr d^ qr , . 



where {ipji}™ f=\ are the Grassmann variables (n variables for each determinant in (II. 6p ). 
Denote 

2 2 
Xj,k = S^iiV'fci + V'w^jOj Xj* = ^W^ju - V'w^'ii): i ^ fc, (2-7) 

2 



to = ^^ji^jh j,k = l,..,n. 



1=1 

-2 



Using that {XjkY = for s > 4, j, k = 1, ..,n (since ipj s = = for any j = 1, ..,n, 
s = 1,2), we expand the second exponent under the integral in (12.61) into the series and 
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integrate with respect to the measure (jl.2p . We get then 

1 - ^ + £j<^) 4 n (i + nn^+,«v. 

/ j=l ^ ' r=l q=l 

where \x± is 4-th moment of the common probability law Q of the entries $iWjk of 

( P . Note that 



where is defined in (11.81) . Thus, (12.81) yields 

D 2 l F 2 (A) = D 2 l / p=i niI d V^ ( 2 - 9 ) 



r=l q=l 



where 



-1 = -^E((^) 2 -fe) 2 )-E(^) 2 (2-io) 



2 

i<fc 3=1 

2 / n \ ^ / n n 



{=1 \j=l / l<Z<s<2 \i=l fc=l 



i<fc \j=i / 



2 



Now we use the formulas 



^^-exp{ab 2 } = J exp{— ax 2 — 2abx}dx 1 (2-H) 

7T f 

— exp{abc} = / exp{— auu — abu — acu}d$t,udQu, 
a J 

where b, c are complex numbers or even Grassmann variables (i.e. sums of the products 
of even number of Grassmann variables), and a is a positive number. For the case of even 
Grassmann variables this formulas can be obtained by expanding the exponent into the 
series and integrating of each term. Therefore, (12.101) - (12.111) imply 



lira 



exp < ——01 



/ exp < -- {J2 T l + 2 Z UabUab ) \ 

£ I \q=l l<a<b<2 J J 



In I 2tt 2 

U% V \q=x i^a^z / j (2.12) 

f] exp J ^ iTpipj^jp + Yl {iucdipjc^jd + iu cd ip jd ipjc) \ dQ, 

j=l I p =l l<c<d<2 



where 



Q = ( - 1 UU ) , (2-13) 



7^2 is the space of 2 x 2 hermitian matrices and d Q is given in (11. 7p for n = 2. Write the 
formula 

exp O^ 2 } = v? / exp ^~i K4 i p2 ^ ri exp | 2p ii^j ip jl ip j2 ip jl 'ip :i2 1 dp (2.i4) 



with 



e(z) = <i X '. x > °' (2.15) 
1 —ix, x < 0. v y 



Substituting (T2~T2|) - Ell in and using (Q - O we can integrate in CTl 
over the Grassmann variables. We obtain 

D 2 1 F 2 (A) = Z 2 j dp J rfge-f trQ2 -'^ 2 (det(Q-zA) + ^^y, (2.16) 

where Q is defined in (12. 13[) and 

A=r^ M, * = i=£ffi. (2.17) 



A 2 ; ' 2tt 2 v ^^ 

Recall that we are interested in A = A + £/np sc (A ), where A = diag{A ,A } and 
£ = diag{£i,£2}- Change variables to Tj — iXo/2 — i£j/np sc (\o) — >• T j, J = 1,2 and note 
that we can move the integration with respect to Tj from line Qz = \o/2+£j/np sc (\o) back 
to the real axis. Indeed, consider the contour Cjr, which is the rectangle with vertices at 
(-R, 0), (-R, A /2 + fa/n), (R, A /2 + ^/n) and (R, 0). Since the integrand in flZZEED is 
analytic in {t 3 -}| =1 , the integral with respect to Tj of this function over CjR is equal to 
0. Besides, the integral over the segments of lines Qz = ±R tends to as R — > oo, since 
the integrand in (I2.16P is a polynomial of Tj multiplied by exp{— nr 2 /2}. Thus, setting 
R — > oo, we obtain that the integral with respect to Tj over the line Qz = Ao/2 + £j/n is 
equal to the integral over the real axis. Hence, we obtain in new variables 

t tr ( Q+%2--^tr(Q+%s-| K4 b 2 - £tr- 



D 2 l F 2 (A) = Z 2 I dp j dQe 2 vv 2 ; S^T"^" 1 " 2 " ^ 2 " p.c(ao) 



x ( det (Q-^) + ^-^ 1 =Z 2 jdpjdQexp^--j^tr(Q + ^ (2.18) 



-l^b'-^tr^l^j^Q) (l + 



«2 

2pe(n 4 ) ' " 



ndet(Q - «A /2) 
where Q is again the hermitian (see (12. 13ft ) and 

/i n (Q) = det n (Q - zA /2)e-ttr(Q+^o/2) 2 _ (219) 

Let qi, q 2 be the eigenvalues of Q. Set 



^n = {(Q,p):a< |a-iAo/ 2 l < A/ = 1,2, |p| < logn}, 
n? = {06%:i< |»-iA /2| < A}. 



(2.20) 
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for sufficiently small a and sufficiently big A (note that if |A | > 5, then \qi — z'A /2| > <5 2 /4 
and we can omit the first inequality in ( I2.20p ). Note that the integral in (12. f 8p over the 
domain max \qi\ > A is 0(e~ nA2 / 4 ), A — > oo and the integral over the domain min \qi\ < a 

lr — .L . 2. L — J. . 2. 

is O(e- nloga_1 ), a 0. If a < \q t - i\ /2\ < A and |p| > logn, then according to fl2^j) . 
(12. 5 p and (I2.13p . the corresponding integral is bounded by 

Z 2 J \»n{Q)\dQ J (I + Cp/n) n e-^ p2 dp = O(e- Clos2n ), (2.21) 

n% bl>logn 

and we can write 



D^F 2 (A) = Z 2 y"e~ ltr(Q+ ^ ) ^- |K4|p2 ~^ tr ^^ +2pe(K4)dct " (Q ^ ) 
x// n (Q) (1 + / n (det(Q - zAo/2),p)) dpdQ + 0(e~ cl ^ n ) 



(2.22) 



where 

/„(det(Q - iAo/2),p) = e-**™**- 1 ®-^ (l + -^^\ - 1. (2.23) 
Note that f n is an analytic function of p and entries of Q, and we have on fl n 

l k 

|/ n (det(Q-zA /2),p)|<^A (2.24) 

n 

where k is independent of n. It is easy to check that 



J ._ J e -|M|p 2 +2p£(M)det- 1 (Q~^) rfp = ^ 71 gM def 2 (Q-iAp/2) + O (g- C log 2 nj 
|p|<logra 

and we obtain that |/| > C 3 > on fi„ (see (12T201) ). Thus, (12T221 yields 

^^(A) = ( _ 1)n 2 2?r2 y MQ) exp |-ztr(Q + 2A /2)£/p sc (A ) 
+ K4 det- 2 (g - zA /2)} (l + 4 1) (det(g - iAo/2))) dQ + O( e - clo § 2 "), 



where 



tr— 2 



/« (det (Q - iAo/2)) = e 2 " - 1 



+ j- 1 e 2 - tr ^kF y e - lKilp2+2p£{Ki)det ~ 1{Q - lAo/2) f n (det(Q -iAo/2), p)dp. (2.25) 

|p|<logn 

According to ( 12.241) . we get that fn \det(Q — iA /2)) is analytic in elements of Q on 
and 

|/«(det(g-zA /2))| < log fc n/n, (2.26) 
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where k is independent of n. 

Let us change variables to Q — U*TU, where U is a unitary matrix and T = 
diagjtx, t 2 }- Then dQ of (jl.7p for n = 2 transforms to (ti — t 2 ) 2 dtidt 2 d fi(U), where 
li(U) is the normalized to unity Haar measure on the unitary group U{2) (see e.g. [9], 
Section 3.3). Hence, since functions det(Q — *A /2) and tr (Q + iA /2) 2 are unitary 
invariant, ( I2.25f) implies 

D^F 2 (A) = J dfi(U) J d^dhfliU-tXo^r (2.27) 

U(2) L£xL£ 1=1 

-^t{tsV^) 2 ~tlU'{TV^)U-^-+n i f\{t r ^)-^ 

x e a=1 r '= 1 
x (l + /«(det(T-^)))+0(e-^ 2 "), 

where 

Lf = {t g R : a < \t -i\ /2\ < A}. (2.28) 

The integral over the unitary group U(2) can be computed using the well-known Harish 
Chandra/Itsykson-Zuber formula (see e.g. [9j, Appendix 5) 

Proposition 1. Let A be the normal n x n matrix with distinct eigenvalues {a,i}2 =1 and 
B = diag{b\, . . . , b n }. Then we have 

J J exp{-^tr(A-U*BU) 2 }A 2 (B)f(B)dUdB 

J exp{~tr(a J -b J ) 2 }^f(b 1 ,...,b n )dB, (2.29) 



U(n) 



n n/2 



n 

where f(B) is any symmetric function of {6j}™ =1; dB = Yl dbj and A (A), A(S) are 

3=1 

Vandermonde determinants for the eigenvalues {cti}™ =1 , {&i}™ = i of A and B. 
Hence, we obtain finally from ( 12.27ft 



* r( - 1) *V«« (2-30) 
e K 4 (*i-a /2)-3( t2 -iAo/2)-2 ^ + /? (2)( T )) dfi ^ + (e- clog2ri ), 

a — 6 



where is defined in ( jX2gD and /,1 2) (T) = /^ 1) (det(T - iA„/2)) is an analytic function 
bounded by log fe n/n if U e L^, I = 1, 2. 

This asymptotic integral representation is used in the section 3 to prove the theorem 
for m = 1. 



2.2 Asymptotic integral representation for F^m- 

Set 

2m 

D 2m = n \A° (n) te), ( 2 - 31 ) 
i=i 



where £> {n) (0 is defined in (Ol) . Using ([22]), we obtain from (JUSD (cf. (ESJ); 

^ A „ rr, ~ ., 2m n 

J qr 



{ 2m n — 9m n 

r e e (>j--ff)j,fc^« W - 



D^E { / '<* ^ «=* (2.32) 



iOtu^ fc 2U? ,— — n ,„ . . 2m_ 



- E — E ® jt *u-4>Hl>ii)- E -# E 2m n _ 

r=l g=l 

where {V^ljjJi are the Grassmann variables (n variables for each determinant). As in 
(12 ,7p we denote 

2m 2m 

xtk = ^i^ji^ki + ipkiipji), Xj,k = ^(^ji^ki ~ TpkM, J + k , (2-33) 



i=l Z=l 
2m 



J.J 

«=1 



Using that (xfk) s = for s > 4m, j, k = 1, .., 2m (since = ip^ = for any j = 1, .., n, 
I = l,..,2m), we expand the exponent under the integral in (I2.32p into the series and 
integrate with respect to the measure (jl.2p . We get then similarly to ( 12. 9 p 



-'" " - 2m *2 P 2m n 



D£F 2m (A)=D£ / e- 1 *=i *=* Hl[d^ qr d^ gr , (2.34) 

where K 2p is cumulants of the probability distribution of entries $lwjk, ^Wjk of (11.21) . i.e. 
the coefficients in the expansion 

i(t) := logE{e^ fe } = ^ ^(ii) q + o(t s ), t 0. 

q=0 V' 

The function <j\ in ( I2.34p is the same as in (I2.10p (but with xf,k of (12.33|) and the sums 
from 1 to 2m instead of from 1 to 2), 

1 9 n 

°* = i!E((^) 4 + (^) 4 ) + ^E(^) 4 ( 2 - 35 ) 

j<k j=l 

n n 

h<l2<si<s 2 j=l k=l 
j n n 

+ I Yl Y^jh'^js 1 i J jhi J js2-^2^ki 1 'ipks 1 ^ki 2 ' l P 

and for p > 3 we have 



h¥ zs i,h^s 2 j=l k=l 



2p 

p (2p)\ ^ yyA i> k > 1 v ±J yA i> k > ' 1 (2p)! ^ 



2m [3] , - n 
,(P) 



E E 4,i E ^jh-^jip^Jjip^s+i-^jhp ■ E ^kh-ipup^ki^s+i-^khp' 

h,..,l 2p =l s=0 j=l k=l 
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where cjj are n- independent positive coefficients and I = (h, hp)- Using (12. lip we have 



n e(«4) 
2n 



T-r j„ E ( u aia 2 i'lJ'2^ia 1 V'j{, 1 %a 2 i>jb 2 +«a 1 a 2 b 1 6 2 ^ja 1 ^ ^a 2 ^ib 2 ) 

|| e a^6i,a 2 ^b 2 - - ^2.36) 

n_ 2£^4) ^ (w Cl c 2 d 1 d 2 ^ cl V'i C2 V'jd 1 ^d 2 +«' C ic 2 d 1 d 2 V'i C i^C2'/'jd 1 V'jd 2 ) 

[j e C i< C2 <d!<d 2 - dWdV, 



where 



dPU = ] [ d$tw h i 2SlS2 d%w h i 2SlS2 , 

h<h<si<S2 

( n \-( 2 D / \ -(2m) 2 (2m-l) 2 



h^si,l 2 ^s 2 

As well, (12~TT]) yields for p > 3 



(2.37) 



/■ (2m if] 

cxp {^o-p} = y exp J -|«2„| ^ 



2m l j j 
(l,..,i 2 p=l s=0 



n-p]^ e B«ifi*^^ <2 ' 38) 

i=i I h,-^2 P =i <?=o 



with / = (ii,..., / 2p ) and 

2m If] , v _[£](2m) 2 P 

= I] IT(i»r,. s f/^r,.«, c; =(-?-) " ,p>3. (2.39) 
ii,.,/^=i-=o ^ |K2 p'/ 

Substituting (I2.36P - ( I2.38P and (I2.12p with sums from 1 to 2m instead of from 1 to 2 in 
(12.341) and using ( 12. 2ft - ( 12. 3 p we can integrate over Grassmann variables in (12. 34ft . We 

get 

D^F 2m (A) = Z m J dQ J dVdRdWe-^ tTQ \(v,w,r)<5> n (iQ + A,v,w,r), (2.40) 

where T-L 2m is the space of hermitian 2m x 2m matrices, 

v = {^aia 2 6i6 2 |ai ^ 61, a 2 7^ b 2 , a x ,a 2 MM = 1, -,2m}, (2.41) 
w = {w aia3 t n h\a 1 <a2<bi<b 2 ,ax,aQ,bi,b2 = l,..,27n}, 
r P = {ri, s \h, ..,l 2p = 1, ..,2m, s = 0, .., [p/2]}, 

= (^3,---,^2m), 

and 

{2m 
— |«4|w — 2|k 4 |Wi/j — |^2p|^ P ^p f • (2.42) 
P=3 J 
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dQ, dV, dR and dW are denned in ([H} for n = 2m, f[2T3Tj) and (JOSH, and 



Ml3 



Wl2 
«23 



U13 
«23 

r 3 



Ul,2m-1 V>2,2m-1 U^^m-l 
\ U2-m,l U2m,2 U2m,3 



■Ui,2m-1 M l,2m \ 

^2,2m-l W2,2m 

^3,2m~l Ws^m 

7~2m-l W2m,2m-1 

M2m-l,2m T~2m 



(2.43) 



is obviously hermitian. We denote also 

A = diag{Ai,A 2 ,...,A 2m }, Z„ 



2 2m 



D 



p=2 



(2.44) 



According to (I2.2p - (12. 3p Q(iQ + A,t>,w,r) in (I2.40p is a polynomial of the entries of 
iQ + A and of {^^s^/™}, {wv 2SlS2 /n}, {n u ..,i 2p ,s/n p/2 } with n-independent coefficients 
and degree at most 2m and 

1. the degree of each variable in $>(iQ + A, v, w, r) is at most one; 

2. $(iQ+A, t>, w, r) does not contain terms C(iQ+A)wi 1 i 2S1S2 /n or C(iQ+A)W/ 1 i 2SlS2 /ra, 

_ _ _ _ 2m _ 

since the terms '4>ji 1 '4>ji 2 '4 , js 1 ' l l , js2 or i ) jhi ) 3hi ) jsii ) js3 cannot be completed to Yl i^ji^ji 

i=i 

only by terms ip^ipjs] 

3. &(iQ + A, f , w, r) can be written as 

2s(k 4 ) , 



$(zQ + A, u, u>, r) = det(zQ + A) 



77 



+ f n (iQ + A, u/n, w/n, r p /n p ' 2 ), (2.45) 



where f n (iQ + A, v/n,w/n,r/n p / 2 ) contains all terms of Q(iQ + A, v, w, r) which are 
0(n~ 3 / 2 ) as n — > oo and as Q, t>, w, r are fixed, and a[ contains linear with respect 
to v terms. In view of (12.31) 



(2.46) 



where q s ,i,p,r is (2w — 2) x (2m — 2) minor of the matrix + A without rows with 
numbers s and / and columns with numbers p and r. 

Recall that we are interested in A = A + £/np sc (\ ), where A = diag{Ao, . . . , Ao} and 
£ = diag{£i , . . ., & m }- Shift now tj - iXo/2 - i£j/np sc (\ ) -)■ r h j = l,..,2m. Then 
similarly to (12. 18ft we obtain in new variables 



D^F 2m (A) = Z m J dQ J 



u n (v, w, r)$ n (iQ + A /2, v, w, r) 



(2.47) 



Hz 



x e 



itr { Q + ^f-iti (q+i^o )f /psc(Ao) _ i tr _4 



*• -Js^^dV dRdW, 



12 



where Q is the hermitian matrix of ( I2.43P and dQ, dV, dR and dW are defined in (jl.7p 
for n = 2m, f l23Tj) and f l2T39|) . The (1) condition of $ yields 



|$(iQ + f , u, w, r)| < Hi 1 + + Ao/2) g .|) II i} + C 



x II ^1 + C 

ai<H2<&l<&2 



n 



2m 2m [2! 



n n n(i+c 

with n-independent C, and (I2.45P yields 

\${iQ + A ,v,w,r)\ < |det(«Q + A /2)| JJ (l + C(Q) 



n,. 



(2.48) 



x 



n { i + c (q) 



2m 2m [fl 



p/2 



n 



p/2 



(2.49) 



n n n(i+c«) 

il<Z2<Sl<«2 p=31i,..,l 2p =l s=0 

Here C(Q) is bounded if a < |g/ — zAo/2| < A, I = 1, .., 2m and {qi}j=i are the eigenvalues 
of Q. Note that if |Aq| > 5 > 0, then |gj — iAo/2| > S 2 everywhere. Denote 

= {{Q,v, w,r) : a < \q s - iX /2\ < A, \v hhsiS2 \ < logn, (2.50) 

|u>Jifa«i<al ^ lo S n > l r Ml ^ lo g™}> 
9% = {Q G "H 2m :a<\q s - i\ /2\ <A,s = l,.., 2m}. 

According to (I2.48P the integral in (I2.47P over the domain max \qi — i\o/2\ > A 

1=1, ..,2m 

IS — > oo and the integral over the domain min \qi — i\o/2\ < a is 

£=l,..,2m 

O(e~ nloga ), a 0. Moreover, the bound (12.491) implies that this integral over the do- 
main, where the absolute value of at least one of {vi^s^}, {wi^s^} or {r ljS } is greater 
then logn but a < \q s — i\o/2\ < A, s — 1, ..,2m, can be bounded by e~ clog n (similarly 
to (12^21"]) ). Therefore, using ([23]) . f l23T|) . and fl2Tl4"|) to bound the integral with \fi n (Q)\, 
we can write 



D 2m F 2m(A) = Z m / /J, n (Q)e 



" tr (Q+^ li )^ T -2 £ ( K 4)det-i(^+Ao/2)<-^tr —g-j, 



(2.51) 



x£ n (v,w,r) (l + / n (Q,t;,u;,r))rfQdKrf J RdH/ + 0(e- clos2ri ), 
where u n , <r^ and v n (v,w,r) are defined in ( 12. 19ft . ( 12.461) and ( 12.421) respectively, and 

fn(iQ + ^f,v,w,r) - 2e(K 4 )a'Jn 



f n (Q,v,w,r) = e^M^W+^K [ i + 



det(zQ + A /2) 



1. (2.52) 



Note that /„ is an analytic function of the entries of Q, and in view of (I2.49P we have on 

\f n (Q,v,w,r)\ Kn-Wloifn, (2.53) 
where k is independent of n. It is easy to check that 



i , n (Q,v,w,r)dV dRdW 



(2.54) 



2m 



J~J^'^-l e K40-(jQ+A /2)det- 2 (jQ+Ao/2) _|_ £^ e -clog 2 n^ 
p=2 
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where 

u n (Q,v,w } r) = exp{-2e(/t4)det _1 (iQ + A /2)a' 1 }u n (v, w, r), 

a{iQ + A /2) = Qh,si,h,s2Qh,s2,h,si ( 2 ' 55 ) 

with qi l!S1 ,i 2 ,s 2 defined in (12.461) (but for the matrix iQ + A /2 instead of iQ + A ). Note 
also that according to the Cauchy-Binet formula (see [6]), we have that a{iQ + A /2) is 
the sum S2 m -2{A) of principal minors of order (2m — 2) x (2m — 2) of the matrix 

A = {iQ* + A /2)(tQ + A /2) = U*(iT + A /2) 2 U, 

where U is a unitary 2m x 2m matrix diagonalizing Q and T = diagjgi, .., q2m}, i-e. 
Q = U*TqU . Since S , 2 m _2(^4) is a coefficient under A 2 in the characteristic polynomial 
det(/l — XI), S2m-2{A) is unitary invariant, and thus a(iQ + Aq/2) is unitary invariant 
too. Therefore, we have on Q n of (12.501) 



a(iQ + ^)det -\iQ + ^) 



1 



and hence |J| > C > 0. This, (I23T|) and f[2~54l) yield 



< a 



D^F 2m (A) 



(2.56) 



x /i n (g) (1 + /W(Q)) rfg + o(e- c1 ^ 2 "), 
where /z„ is defined in (I2.19p and 
-itr— ? 



f<P(Q) = e »>' 



-^tr— £ 



+ J- i e 2n ^ / u n (Q,v,w,r)f n (Q,v,w 1 r)dVdRdW (2.57) 
n n 



with J of (12.541) and z/ n of ( 12.55(1 . According to ( 12.52(1 and bound from below of |/| on £l n , 
we have 

\f£\Q)\< Vn/n 1 / 2 . (2.58) 

Besides, fn\Q) is analytic in elements of Q. 

Let us change variables to Q — U*TU, where U is a unitary 2m x 2m matrix and 
T = diagjti, . . . , t^m}- The differential dQ in (I2.56P transforms to A 2 (T)dTd n(U), where 

2m 

dT = Y\ dti, A(T) is a Vandermonde determinant of {t;} 2 ^, and n(U) is the normalized 



i=i 



to unity Haar measure on the unitary group U(2m) (see e.g. j9], Section 3.3). Functions 
det(iQ + 4p), tr(Q + ^-) 2 and cr(iQ + 4f) (as we proved before) are unitary invariant. 
Hence, ( 12. 56(1 implies 

/ i \mn„2m 2 r 



2m 



dTjliU 



iA , 



U(2m) 



(^) 2 



xe 



(2.59) 



l<Ks<2m 



x A 2 (T) (l + fP(U*TU)\ + 0(e- cl °s 2 
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where is defined in (I2.28P . Using Proposition [TJ we have 



(_l)mn n 2 m 2 r _» g (ti + ^fl)2 +K4 £ ( ti _^a )-2 (t ._iM r 2 

1 / p ( = 1 Ki<j<2ra 



2m 



D 2m 2 m n 2m2 

{Li) 2 

( -/fft^)^- \ (2 ' 60) 



\2m 



A 2 (T)n(t ( -f 



+ f (2) (T) 
v A(T)A(-ze/p sc (A )) n j 



dT + 0(e 



-clog n\ 



where 



ff\T) = J exp|-.tr^(T + ^)[/-^|/( 1 )(f/W)rfM^)- 
According to ( 12. 58ft . we get that 

\J^\T)\ < n-^log'n^ eLi,l = 1, ..2m. (2.61) 
Hence, we obtain finally 

-?Eft4) 2 -E(t ! 4)^+^ E (t h -^)'Ht l2 -^)- 2 (2.62) 

l=i l=i h< l 2 

2m 



x (l + /( 2 )(T)) [I dtj + Oie 



-clog 2 n"i 



where 

2m 

/f (T) = A(T)A(-^„(Ao)) e , S' ,S ' / '" (A - ) /f (IT). 
fn\T) is an analytic function bounded by tit 1 / 2 log fc n if G L^, / = !,.., 2m. 



3 Asymptotic analysis. 

In this section we prove Theorem [1] passing to the limit n — )■ oo in f |2.62j) for Aj = 
A + £j/np sc (\ ), where p sc is defined in ( Oil . A G (-2,2) and G [— M, M] C E, 
j = l,..,2m. 

To this end consider the function 

V(t, A ) = t + l -h t - log(f - zAo/2) - i^i. (3.1) 

Then (127621 yields 

r-^^A)^™, / H/ ri (t 1 ,...,t 2m )rfT + 0(e- clos2n ), (3.2) 

w sc (A )) m y 

(L A )2m 
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where Z) 2m is defined in (I2.3ip . 



A(0 (3.3) 



xe l<Ks<2m " ~ (l+fP{T) 



and 



/ -\\mn ri m 2 n (\ \m(m-l)„-mK4 

m > n ~~ (•_^m(2m-l)22m 7r 2m 2 ' ^ > 



Now we need 

Lemma 1. The junction ?R.V(t, Aq) for teK /ias t/ie minimum at the points 



t = X± :=±^5. (3.5) 

Moreover, ift ^ U n (x±) := n -1 / 2 logn, x ± +n~ 1//2 logn), then we have for sufficiently 
big n 

nV(t,X)> Cl °f n . (3.6) 

Proof. Note that for t E R 

W(t, A ) = i (t 2 - (4 - A 2 )/4 - log(t 2 + A 2 /4)) , (3.7) 



and thus 

d 2 1 2t 



-^(t,A ) = t- - ■ x2//| , (3.f 



r^(t,A ) = l- 7^o77 + 



dt 2 v ' uy t 2 + A 2 /4 (t 2 + A 2 /4) 2 ' 

Therefore, t = x± of (13. 5p are the minimum points of ?R.V(t, Aq). Note that 



V + :=V{x + ,Xo) = ZA °^ A ° -zarcsin(-A /2), (3.9) 



^_:=y(x_,A ) = - ZA °^ A ° + zarcsin(-A /2) - m. 

Thus we have 

&V(x±,\ ) = 0. 

Expanding 3lV(t, A ) into the Taylor series, we obtain for t G U n (x±), using (13. 8p - (13. 9 p 

^(t,A ) = ^^(t-x±) 2 + 0(n- 3/2 log 3 n), (3.10) 
where x± is defined in (13 .5h . Hence, for t G" t/„(a;±) we get 

KV(t,A)>^^, 
n 

which completes the proof of the lemma. □ 
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Next note that since \tj — i\o/2\ > a for tj G L^, j = 1, .., 2m, we have 



cxp 



L 4 J2 (t l -i\ /2y 2 (t s -i\ /2)- 2 \ 

I Kks<2m J 



< C. 



(3.11) 



This, the inequality |A(T)/A(f)| < G x for \tj\ < A, j = 1, . .,2m and distinct {£i)j=i, 
(EED and fl3TTTD yield 



Z,,,,, / / / W„(*i, . . . ,t 2m )dT < C in m V C2log2n , 

LA\(U+\JU.)LA L A 

where , W n and Z mj „ are defined in (I2.28p . (13. 3 p and (I3.4p respectively, and 

U± = {teR:\t-x±\<n- 1/2 logn} (3.12) 



with x± of (I3.5p . 

Note that we have for t G £/± in view of (13.11) and (13. 9p as n -> oo 



V(t,A ) = V±+ 1 



1 \ (t-x±) J 



(x ± - i\ /2f 



f±(t-x ± ), (3.13) 



where f±(t — x±) = 0((t — x±) 3 ). Shifting tj — x± — > tj for tj G U± we obtain using ( 13. 9ft 
that the r.h.s. of (13.21) can be rewritten as 



2m 2m 



/ n^ii e 2 psc(Ao) /qj(j) — 



an • • ■ i ^2m "I - -£a2 



2m J 



2m 

«4 E (^+Pai)~ 2 (*s+Pa s )" 2 - E (nV'a ■ +ix Q £j/p S c(A )) 
3 KKs<2m J = l 



(3.14) 



(l + /i 2) (T)) + 0(e- cl °^"), 
where sum is over all collection a = {ctj} 2 ^, atj = ±, j = 1, .., 2m and 

c± = l+p± 2 , p± = x±-i\ /2, U n = (-2n- 1/2 logn,2n- 1/2 logn). 
Note that 

2m 



(3.15) 



/ := 



A(ti + x ai , . . . ,t 2m + x a2m )Y[dtj 



3=1 
2m 



det 



"P sc (A )c o 



fc-i 



(3.16) 



dtj 



j,k=l 



where 



U hn = \-2n- l l 2 \ogn + 



Kj n_-l/2i , 



^P S c(Ao)c a? 



2n 1/2 log n + 



np sc (X )c o 
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Since f±(t) = 0(t 3 ), changing variables to y/ntj — > tj, expanding ex^>{—nf a .{tj/^/n — 
i^j/np sc (\o)c a )} in (I3.16p . and keeping the terms up to the order n -4 ™ 2 , we obtain as 
n — t- oo 

iJ uM det { ( x " - dfef + ^ Wb, C (i + <3 - i7) 

where Pj^l and Pf^l are polynomials with n- and j-independent (but fc-dependent) coef- 
ficients of degree at most 4m 2 . Consider 

{fc-i 1 2m 

K ~ - m V ) + AP M^» " ( 3 - 18 ) 

Note that _D(£/n, A) is a polynomial of {£j/ n Psc(Ao)} 2 =i and A. Without loss of generality, 
let «i = ... = a s = +, = ... = a 2m = — • Then it is easy to see that if £ 3 - = for 
j, / = 1, .., s or j, I — s + 1, .., 2m, then D(£/n, A) = 0. Thus, 

£>(£/n, A) = Afo/n, . . . , £» A(£ s+1 /n, . . . , 6m/n)(C + AF(£/n, A)), (3.19) 

where F(£/n,\) is a polynomial with bonded coefficients. Substituting A = in (I3.18P 
and computing the Vandermonde determinant, we obtain 

s(s — 1) (2m-s)(2m- 3 — 1) 

2 / o \ 2 

C = " 

2m 



p sc (A )c + y VP^(A )c 



if* 



n n [x + -x-- — ^ + 

j=l fc=s+l 



np sc (A )c+ np sc (A )c 

s(s — 1) (2m-a)(2m — 3-1) 

— i 



(x + -x_) s ( 2 — s )(l + o(l)). 

x p sc (A )c + y vPs C (A )c_y 

Hence, for ai = ... = a s = +, = ... = a2 m = — we get from (I3.19P as n — > oo 

n m2 J _ 2 m vr m (-i/p sc (A )) ?n(m ~ 1)+{m - s)2 n-( m - s ) 2 ( ^4 - Ag) s ( 2m " s )(l + o(l)) 

s 2^ • ( 3 ' 2 °) 



j=lZ=s+l 



(c+)sV 2( C _)(2 m - s) v2 n n 



This expression has the order at most 0(1), and for s 7^ m it is of order o(l). Hence, the 
terms of ( 13. 14ft are always of order 0(1) and the equality holds only if m of {ctj} 2 ^ are 
pluses, and m last ones are minuses. Consider one of such terms in ( I3.14p . for example 
«i = ..= a m = 1, a m+ i = .. = aim — — 1- Substituting the expressions (13. 4p . (13 . 1 5[) and 
( I3.20p with s = m we can rewrite this term as 



a(m-l)K 4 (A§-2) 2 /2 ^OT(m+l)giir($m+l+-+&m— Si — ••— ^m) 

7r 2m 2 -2m m 
(JWtt)™ EI 
»J=1 



(3.21) 



In view of identity 



A(f 1 ,..,£ m )A(f m+1 ,..,f 2m ) (2z7r) m A(ei,..,e m )A(e m+ i,..,6 m ) 
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2m 



the determinant in the l.h.s. is the sum of exp{i7r ^2 otj^j} over the collection {ctj}^, in 

which m of elements are pluses, and m last ones are minuses, with certain coefficient. In 
view of the identity (see [12], Problem 7.3) 



m(m-l) k< l 
•1) 2 



n(a fc - ai)(b k - k 



det 



II ( a k ~ k) 



1 



3 J fcj=l 



fcj=i 



the coefficient under exp{i7r(£ m+1 + .. + £ 2m — £i 

1 



det 



J ) j,k=l 



Cm)} is equal to 



, 1 x m(m-l) 

— 1 2 



(2i7r) w A(a, .., £ m )A(W, ■■, 6m) 



»J=1 



Other coefficients can be computed analogously. Thus, restricting the sum in (13.14p to 
that over the collection {otj}*™^ in which m of elements are pluses, and m last ones are 
minuses, and using (I3.2ip . we obtain Theorem [1] after certain algebra. 
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